Extracting neutron generalized parton distributions from 3 He data 
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An impulse approximation analysis is described of the generalized parton distributions (GPDs) H 
and E of the 3 He nucleus, quantities which are accessible in hard exclusive processes, such as coherent 
deeply virtual Compton scattering (DVCS). The calculation is based on the Avl8 interaction. The 
electromagnetic form factors are correctly recovered in the proper limits. The sum of the GPDs H 
and E of 3 He, at low momentum transfer, is largely dominated by the neutron contribution, thanks 
to the unique spin structure of 3 He. This nucleus is therefore very promising for the extraction of the 
neutron information. By increasing the momentum transfer, however, this conclusion is somehow 
hindered by the the fast growing proton contribution. Besides, even when the neutron contribution 
to the GPDs of 3 He is largely dominating, the procedure of extracting the neutron GPDs from 
it could be, in principle, nontrivial. A technique is therefore proposed, independent on both the 
nuclear potential and the nucleon model used in the calculation, able to take into account the nuclear 
effects included in the I A analysis and to safely extract the neutron information, even at moderate 
values of the momentum transfer. Thanks to this observation, coherent DVCS should be considered 
a key experiment to access the neutron GPDs and, in turn, the orbital angular momentum of the 
partons in the neutron. 
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INTRODUCTION 

Generalized Parton Distributions (GPDs) (Tj-Q 
parametrize the non-perturbative hadron structure in 
hard exclusive processes, encoding therefore a wealth of 
information. For example, the hadron three-dimensional 
structure [4( and the parton total angular momentum 
could be unveiled by the measurement of GPDs, which 
will be therefore a major achievement for Hadronic 
Physics in the next few years. In particular, in order 
to access the hadron angular momentum content, from 
which the orbital angular momentum (OAM) part 
could be estimated by subtracting the helicity one, 
measurable in deep inelastic scattering inclusive (DIS) 
and semi-inclusive (SiDIS) processes, the knowledge of 
two GPDs is mandatory [3]. They are the two chiral 
even, parton helicity independent GPDs occurring 
at leading twist, i.e., H, a target helicity-conserving 
quantity, and E, a target helicity- flip one. The cleanest 
experiment to access them is Deeply Virtual Compton 
Scattering (DVCS), i.e. the process eH — > e'H'^f when 
Q 2 ^> M 2 (everywhere in this paper, Q 2 = — q ■ q is the 
momentum transfer between the leptons e and e', A 2 
the one between the hadrons H and H' and M is the 
nucleon mass) 0, Hj]. Despite severe difficulties related 
to the complicated way GPDs enter the measured cross 
sections, DVCS data for proton and nuclear targets axe 
being analyzed (recent results can be found in Refs. 
and GPDs are being extracted from them (see Refs. 



and references therein). 

The measurement of GPDs for nuclear targets has been 
proposed in several papers, for the deuteron [8|, |9[ and 
for finite nuclei 0, |To|. One of the main observations 
has been that, knowing GPDs, one could unveil medium 
modifications of bound nucleons, otherwise difficult to 
distinguish from other medium effects in data collected 
through standard DIS experiments. These investigations 
could therefore seriously contribute to shed some light on 
the origin of the so called EMC effect [11], a puzzle still 
far to be solved. 

Great attention has anyway to be paid to avoid con- 
fusing unusual effects with conventional ones. To this re- 
spect, few-body nuclear targets, for which realistic stud- 
ies are possible and exotic effects are in principle distin- 
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guishable, play a special role. To this aim, in Ref. [121 ] . 
an Impulse Approximation (I A) calculation of i7 3 , the 
GPD of 3 He corresponding to the flavor q, has been pre- 
sented, valid for A 2 <Q 2 , M 2 . The approach permits to 
investigate the coherent, no break-up channel of DVCS 
off 3 He, which cannot be observed at large A 2 , due to 
the vanishing cross section. The main conclusion was 
that the nuclear GPDs cannot be trivially inferred from 
those of nuclear parton distributions (PDFs), measured 
in DIS. 

In a recent Rapid Communication of ours, Ref. [l3| . 
the approach of Ref. [12j has been extended to evaluate 
the GPD E q of 3 He, E^. The main goal was to study the 
possibility of accessing the neutron information, which is 
very relevant because it permits, together with the proton 
one, a flavor decomposition of GPDs data. One should 
not forget that the properties of the free neutron have 
to be investigated using nuclear targets, taking nuclear 
effects properly into account. In particular 3 He, thanks 



2 



tron, such as its helicity-flip GPD E q 
been found that the sum of the GPDs Hj* 



to its peculiar spin structure, has the unique property 
of simulating an effective polarized free neutron target 
(see, i.e., [1 



16|). 3 He is therefore a serious candidate 
to measure the polarization properties of the free neu- 

In Ref. [13] it has 
md E^, at 

low momentum transfer, is indeed dominated to a large 
extent by the neutron contribution, making 3 He targets 
very promising for the extraction of the neutron infor- 
mation. Unfortunately this is not the end of the story. 
The same analysis has shown in fact that the proton 
contribution grows fast with increasing the momentum 
transfer. Besides, even if the neutron contribution to the 
GPDs of 3 He were largely dominating, the procedure of 
extracting the neutron GPDs from it could be, in prin- 
ciple, nontrivial. In this paper, a more comprehensive 
analysis is presented. In particular, a technique is pro- 
posed, independent on both the nuclear potential and the 
nucleon model used in the calculation, able to take into 
account the nuclear effects included in the IA analysis 
and to safely disentangle the neutron information from 
them, even at moderate values of the momentum trans- 
fer. Thanks to this observation, coherent DVCS off 3 He 
is strongly confirmed as a key experiment to access the 
neutron GPDs. 

The paper is structured as follows. In the next section, 
the formalism used in the I A analysis, only sketched in 
the previous Rapid Communication, Ref. [13], is devel- 
oped and motivated. In the third one, the ingredients 
used in the calculation are described. In the following 
section, the proton and neutron contributions to the 3 He 
GPDs, obtained in the present approach, are presented, 
together with their integrals, providing correctly the elec- 
tromagnetic form factors (ffs). In the fifth section, it 
comes the most relevant result of the paper: a safe ex- 
traction procedure of the neutron information from 3 He 
data, independent on both the nuclear potential and the 
nucleon model used in the calculation. The impact of this 
study in the present experimental scenario is eventually 
discussed in the Conclusions. Two Appendixes have been 
added with some details on the description of 3 He GPDs 
in terms of few-body wave functions. 



through the light cone correlator 
F^(x,A 2 ,0 = J ^e^ p+ ^ A {P's'\d^ q \Ps) A 



Z +=0,Z±=0 



2P+ 



Hq(x, A 2 , £)^(P', s')^u{P, s) 



Anna. A 

E^x,A 2 ,Ou(P',s') t -^^u(P,s) 



where C>£ = % (-§) ^ 4> q (§), P = (P+P')/2, % is the 
quark field, m is the hadron mass and = (qo,q). El- 
lipses denote higher twist structures. The £ variable, the 
so-called skewedness, is £ = — A + /(2P + ) (everywhere in 
this paper, a ± = (a ± a 3 )/^)- In addition to the vari- 
ables x, £ and A 2 , GPDs depend on the momentum scale 
Q 2 . Such a dependence, irrelevant in this investigation, 
is not shown in the following. 

In proper limits, GPDs are related to known quantities, 
such as the PDFs and the electromagnetic form factors 
(ffs). In particular, the following constraints, relevant for 
what follows, hold: 

i) in the "forward" limit, P' = P, where A 2 = £ = 0, 
DIS Physics is recovered, and H^(x, A 2 , £) coincides with 
the usual PDF, Hf(x,0,0) = q(x), while E£(x,0,0) is 
not accessible] 

ii) the integration over x gives, for H q (E q ), the con- 
tribution of the quark of flavour q to the Dirac (Pauli) ff 
of the target A: 



dxHf(E?)(x,tf,Z) = Ff«(A 



(2) 



iii) the polynomiality property, involving higher mo- 
ments of GPDs, according to which the x-integrals of 
x n H q and of x n E q are polynomials in £ of order n + 1. 

For later convenience, let us define the following auxil- 
iary function, given simply by the sum of the GPDs 
and E^ for a given target A of spin 1/2: 



Gii"(x,A 2 ,0 = H q A (x,A 2 ,0 + E q A (x,A 2 ,0 . (3) 

This function, due to Eq. ((2]), fulfills obviously the 
following relation 



(1) 



IMPULSE APPROXIMATION ANALYSIS OF 3 HE 
GPDS 

Let us first remind the main properties of GPDS (for 
comprehensive reviews, see, e.g., [17H19|). to really un- 
derstand the importance of measuring the neutron GPDs 
and the advantages offered by 3 He. 

For a spin 1/2 hadron target, with initial (final) 
momentum and helicity P(P f ) and s(s'), respectively, 
the GPDs i^(x,A 2 ,£) and E q (x, A 2 , £) are introduced 



J dxGfi q (x,A 2 ,Z) = F^\A 2 ) + F^{A 2 ) 

= Gii"(A 2 ) , (4) 

being G^ q (A 2 ) the contribution of the quark of flavour 
q to the magnetic ff of the target A. 

A fundamental result is Ji's sum rule (JSR) [3], accord- 
ing to which the forward limit of the second moment of 
the unpolarized GPDs is related to the total angular mo- 
mentum of the quark q in the target A, (J q ), according 



3 



to 



J dxxG^ q (x,0,0) 



(5) 



The combination G^J q 



Eq is therefore needed 



to study the angular momentum content of the nucleon 
iV, through the JSR, and OAM could be obtained from 
(Jq), being the helicity content measurable in DIS and 
SiDIS. The measurement of GPDs is therefore very help- 
ful for the understanding of the cumbersome spin struc- 
ture of the nucleon. The proton data alone do not allow a 
clear flavor decomposition of the GPDs and therefore, as 
it happens for any other parton observable, the neutron 
data are very important. To obtain the latter, 3 He is, 
in principle, an ideal target. In facts this system, among 
the light nuclei, is the only one for which the function 
G^ q could be dominated by the one of the neutron, be- 
cause: i) 4 He is scalar and it does not show up any 
ii) 2 H is useless to this aim. To explain the latter state- 
ment, let us think to the forward limit, relevant for the 
JSR. According to Eq. (2) the size of for a given 



target A can be related to its ffs, whose normalization 
is 2^(0) = Z A , Ff(0) + F 2 A (0) = G&(0) = ha, being 
Za and \±a the charge and the magnetic dipole moment 
of the target, respectively. Using the experimental data, 
one gets, for 2 H, F| (0) ~ — 0.14/i/v, reflecting a small E 2 q 
(in the analysis of Ref. [8[ , this contribution has been in- 
deed neglected). On the contrary, in the 3 He case, F 3 (0) 
is sizable {p± — 4.13/i/v), and, moreover, if summed to 
F-f (0), yields \±<$ ~ — 2.13/iiv, a value rather close to the 
neutron one, fi n ~ — 1.91/i/v- As it is well known, [13 
and fjL n would be equal, i.e., there would be no proton 
contribution to /i3, if 3 He could be described by an in- 
dependent particle model with central forces only. Of 
course this scenario is a rough approximation; nonethe- 
less, realistic calculations show that the system lies in this 
configuration with a probability close to 90 %, allowing 
a safe extraction of the neutron DIS structure functions 
from 3 He data, as suggested in (TH, [l6|, estimating effec- 
tively nuclear corrections using static properties. Here, 
the situation is in principle somehow different, because 
GPDs are not densities. Anyway, this scenario is recov- 
ered at least in the forward limit, where the JSR holds: 
static 3 He properties can be again advocated. 

In a previous Rapid Communication of ours, Ref. [l3| . 
it has been established to what extent, close to the for- 
ward limit and slightly beyond it, the measured GPDs 
of 3 He are dominated by the neutron ones. Ref. [Hj 
represents a pre-requisite for any experiment of coherent 
DVCS off 3 He, an issue which is under consideration at 
JLab. In this section the formalism used in the I A analy- 
sis, only sketched in the previous Rapid Communication, 
Ref. [13], is developed and motivated. 

The GPD Hq of 3 He has been obtained, in IA, already 
in Ref. 12L Let us see now how the scheme can be gen- 

~3,q 
M J 



Eq. ([3]). First of all, one should notice that, in addition 
to the kinematical variables x and £, one needs the cor- 
responding ones for the nucleons in the target nuclei, x' 
and The latter quantities can be obtained defining 
the "+" components of the momentum k and k + A of 
the struck parton before and after the interaction, with 
respect to p + — \(p + p') + , being p(p f ) the initial (final) 
momentum of the interacting bound nucleon [12j: 



(k- 



k+ 

A)+ 



(x - 
(x 



(x' 
(x' 



e)p + , 
e)p + , 



so that 



A+ 

2p+ 

X = —X , 



and, since £ = -A+/(2P+), if z = p+/P+ 



2(i + -£ 



(6) 
(7) 



(8) 
(9) 

(10) 



Now, the standard procedure developed in IA studies 
of DIS off nuclei (see, i.e., [20]) is applied to obtain, for 
i7 3 , G^?, convolution-like equations in terms of the cor- 



qi w M 

responding nucleon quantities, , G^J q . Let us just 
recall the main steps of the derivation. First of all, since 
the nuclear states have to be described by non-relativistic 
(NR) wave functions, the nuclear overlaps and the one- 
body operator in Eq. (pp) are treated in a NR manner; 
then, proper components of the 7 M matrix and proper 
combinations of the nuclear and nucleon spin projections 
are selected to extract, from the NR reduction of the cor- 
relator Eq. (pQ), independent relations for i^ 3 , G^?. This 
procedure is detailed in Appendix A. The final result is, 
for any spin 1/2 target A (cf. Eq. JTJ): 

H?(x,A\0 = ^Flt°(x,A*,0, 
1 m 



2P + 



F^\x,A\i). (11) 



eralized to obtain also the combination of GPDs G 3,q 



Once the correct NR treatment is identified, IA is ap- 
plied. The detailed machinery is thoroughly explained in 
Ref. [12] for the case of and it is not repeated here. 
The result, already presented in Ref. [13j, is ob- 
tained analogously. The procedure requires the insertion 
of complete sets of states to the left and to the right hand 
side of the one-body operator O q in Eq. (pQ), so that one- 
body matrix elements and nuclear overlaps are identified 
thanks to the IA. Using then Eq. (fTT|) to obtain rela- 
tions between nuclear and nucleon GPDs, from the ones 
between nuclear and nucleon correlators, the following 
convolution-like formulae are eventually found: 



4 



H*(x,A 2 ,0 = J2f dE f dpY.Y. P ss,s.Mp',E) 

AT J J Co 

(12) 



S s 



N ^ 



dp 



x [P»_ t+ _tfy, E) - P»_ t _ + {p,p', E)} 



(13) 



In Eqs. ([T2]) and (fT3]h proper components appear of the 
spin-dependent non-diagonal spectral function of the nu- 
cleon N in 3 He: 



1 mVme 



Pg s >,ss>(P,P , ,E) = (2 ; )6 - V 2 - / dSl t 



(14) 



Y,{P'S'W's',ts t ) N {ps,ts t \PS) 



N 



where S,S'(s,s') are the nuclear (nucleon) spin projec- 
tions in the initial (final) state, respectively, and E = 
E min + E^, being E* R the excitation energy of the two- 
body recoiling system and E min = \E3 He \ — \E2 H \ = 5.5 
MeV. The most important quantity appearing in the def- 
inition Eq. ([T4]) is the intrinsic overlap integral 



(ps,ts t \PS) N = / dye^ix 



Cat? 



^(£)|*f(z,$) 



(15) 



between the wave function of 3 He, , with the final 
state, described by two wave functions. One of them is 
the eigenfunction \f^% with eigenvalue E = E min + E R , 
of the state s t of the intrinsic Hamiltonian pertaining 
to the system of two interacting nucleons with relative 
momentum t, which can be either a bound or a scattering 
state. The other one is the plane wave representing the 
nucleon N in IA. The relation between the overlap Eq. 
([T5]) and the 2- and 3-body radial wave functions at our 
disposal is discussed and explained in the Appendix B. 
As discussed in Ref. [12j, the accuracy of this calcu- 



lation, since a NR spectral function is used to evaluate 
Eqs. (H3 and is of order O (p 2 /M 2 , A 2 /M 2 ) . The 

interest of the present calculation is precisely that of in- 
vestigating nuclear effects at low values of A 2 , for which 
measurements in the coherent channel may be performed. 

In order to evaluate Eq. (fl3]h the quantity of inter- 
est here, one needs models of the nuclear overlaps, Eq. 
(USD, and of the nucleon GPDs and . The models 
used in the present calculation are described in the next 
section. 
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FIG. 1: (Color online) The magnetic ff of 3 He, G 3 M (A 2 ). Full 
line: the present I A calculation, obtained as the integral of 
Eq. (|13[h summed over the flavors (see text), together with 
the experimental data [32]. The units used for the momen- 
tum transfer, A M = V— A 2 , have been chosen for an easy 
comparison with the results of Ref. [30] • 
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FIG. 2: (Color online) As in Fig. 1, but in units of GeV 2 , as 
a function of —A 2 . 



INGREDIENTS OF THE CALCULATION 

In order to evaluate Eq. (fT3|) . the exact nuclear over- 
laps described above, Eq. (fT5]h have been obtained along 
the line of Ref. [21[ , using the wave function of Ref. [22| 
corresponding to the Avl8 interaction [23|, taking into 
account the Coulomb repulsion between the two protons. 
For a relevant comparison to be described later on, the 
calculation has been performed also using the nuclear 
overlaps corresponding to the Avl4 interaction [24[, in- 
cluding again the Coulomb repulsion, firstly evaluated to 
be used in Ref. [21] . 

For the nucleonic GPDs, the model of Ref. [25| has 
been used, which, despite of its simplicity, fulfills the 
general properties of GPDs. The model for the GPD 
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H q is built in agreement with the Double Distribution 
representation [26|: 



-1 



#f(x,A 2 ,£) = 



5(5 + £a — x) 



x 4> g (x, A 2 ,a)da , (16) 

using the factorized ansatz: 

4> g (x, A 2 , a) = h q (x, a)$ q (x)F q (L\ 2 ) . (17) 

The expressions for the functions h q (x, a) and <& q (x) can 
be found in [l2j; we recall here that the F q (A 2 ) term in 
Eq. (fT7]h i.e. the contribution of the quark of flavour g to 
the nucleon ff, has been obtained from the experimental 
values of the proton, Ff , and of the neutron, Ff, Dirac 
ffs. 

For the u and d flavors, neglecting the effect of the 
strange quarks, one has 



i[2Ff(A 2 )+F»(A 2 )] 



F U (A 2 ) 

F d (A 2 ) = 2F 1 "(A 2 )+Ff(A 2 ) . 



(18) 



The contributions of the flavors u and d to the proton 
and neutron ff are therefore 



^(A 2 ) 



r d 



f u (A 2 ) , 



(19) 
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FIG. 3: (Color online) The quantity x^& M {x^ A 2 , £), shown 
in the forward limit (full), together with the neutron (dashed) 
and the proton (dot-dashed) contribution (after Ref. [131]). 



THE NEUTRON AND PROTON 
CONTRIBUTIONS TO THE GPDS OF 3 HE 

Using the ingredients described above, the calculation 
of Eq. ([T3]) has been performed in the nuclear Breit- 
frame. Unfortunately, the only safe way to establish the 
reliability of the approach, the comparison with experi- 
ments, is not possible here. In facts, data for the GPDs 
are not available and, for in particular, even the for- 
ward limit is unknown. One check is in any case possible 
and it is therefore a crucial one: the quantity 



and 



(21) 



F U "(A 2 ) 
F d "(A 2 ) 



^(A 2 ) , 
-^«(A 2 ) , 



(20) 



respectively. 

For the numerical evaluation, use has been made of the 
parametrization of the nucleon Dirac ff given in Ref. [I?} . 

The model has been minimally extended to 
parametrize also the GPD E q , assuming that it is 
proportional to the charge of q (this choice, a very 
natural one, is used, e.g., in Refs. [28]). In this way, the 
same relations ([T9]) , ([20]) are found between the u and d 
contributions to the Pauli ff and the proton and neutron 
Pauli ff; for the latter quantities, the parametrization of 
Ref. [27j has been used. 

For a relevant comparison to be described later on, the 
calculation has been performed also using a very different 
model of the GPDs, i.e., the one described in Ref. (29| . 
arising in a constituent quark description of the nucleon 
structure, performing a microscopic calculation without 
assuming any factorized form. That model refers to a low 
Q 2 scale and it is reasonable only in the valence quark 
region. 



i.e., Eq. (|3j), summed over the active flavors, can be in- 
tegrated over x to give the experimentally well-known 
magnetic ff of 3 He, Gf^A 2 ) (cf. Eq. flU). As we have 
already written in Ref. [13], the result obtained using 
this procedure is in excellent agreement with the Avl8 
one-body calculation presented in Ref. [30] , and with the 
non-relativistic part of the calculation in Ref. 31] (see 



Figs. 1-2). Moreover, for the values of A 2 which are rel- 
evant for the coherent process under investigation here, 
i.e., —A 2 < 0.15 GeV 2 , our results compare well also 
with the data [32]. For higher values, the agreement is 
lost and to get a good description one should go beyond 
IA, including three-body forces and two-body currents 
(see, e.g., [30]). If measurements were performed at high 
values of —A 2 , our calculations could be improved by 
allowing for these effects; for the moment being, since 
coherent DVCS cannot be measured at high — A 2 for nu- 
clear targets, the good description obtained close to the 
static point is quite satisfactory. 

With the comfort of this successful check, let us briefly 
recall the main outcome of Ref. [13], having a look at 
the nuclear GPDs. The quantity x^& M {x^ A 2 , £), which, 
in the forward limit, yields the integrand of the JSR 
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FIG. 4: (Color online) The same as in Fig. 3, but at A 2 
-0.1 GeV 2 and 6=0.1 (after Ref. [ljj). 



The quantity x^G^{x^ 



A 2 ,0 for the d (full) and 
-0.1 GeV 2 , and $3=0.1. The 
neutron contributions for the d (dashed) and u (long-dashed) 
flavor are also shown (after Ref. [131]). 



FIG. 6: 

u (dot-dashed) flavor, at A 2 
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FIG. 5: (Color online) The proton contribution, G^ P (A 2 ) 
(dot-dashed), and the neutron one, G 3 ^ n (A 2 ) (dashed), to 
the magnetic ff of 3 He, G^(A 2 ) (full), obtained within the 
present I A approach. 



(cf. Eq. (|5]) where the relation is given for a given fla- 
vor q), is shown in the forward limit (Fig. 3), and at 
A 2 = -0.1 GeV 2 and £ 3 = 0.1 (Fig. 4). Here we have 
defined x 3 (M A /M)x and £ 3 (M A /M)£, in order 
to recover the standard notation used in studies of DIS 
off nuclei. In facts, the Bjorken variable is defined as 
%Bj — Q 2 1 '(2Mv), being v = qo in the laboratory sys- 
tem. It ranges naturally between and Ma/M for a 
nuclear target of mass Ma- It is therefore convenient to 
rescale the variables by the factor Ma/M. 

The difference in size between the curves in Figs. 3 and 
4 reflects the effect of increasing A 2 , a dramatic behavior 
basically governed by that of the ff. This fact can be re- 
alized looking at Fig. 5, where the separate contribution 
of the neutron and the protons to G\ I (A 2 ) is shown in 
linear scale. 

The most striking result is actually that the contri- 
bution of the neutron is impressively dominating the nu- 




-0.2 



FIG. 7: (Color online) The ratio of the proton to neutron 
contribution to the quantity xzG 3 M (x, A 2 , £) of 3 He (full), Eq. 



(T23T) . obtained at -A z = 0.1 GeV and £ 3 = 0, using the 
nucleon GPDs model of Ref. [2^] in Eq. ([T5|h compared with 
the same ratio evaluated using the nucleon GPDs model of 
Ref. [H (dashed). 



clear GPD at low A 2 , with the proton contribution grow- 
ing fast with increasing A 2 . However, as shown in Fig. 
6, for the flavor d the impressive dominance of the neu- 
tron contribution varies slowly with increasing A 2 . These 
facts can be understood qualitatively looking at Eq. (fT3|) 
which, being rather involved, can be usefully sketched as 
follows 



G 



3,q 
M 



(22) 



where P^ n ^ describes the proton (neutron) dynamics in 

3 He, while G P 1 ^ n - >,q is the contribution of the flavor q to the 
GPDs of the free proton (neutron). As already explained, 
due to the spin structure of 3 He, P 3 is quite larger than 
P 3 justifying the relevance of the neutron contributions 
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in Fig. 3. Anyway, with increasing A 2 , for the u flavor, 
the term G^ n gets much larger than G^} n , explaining 
the growth with A 2 of the relative size of the proton 
contribution with respect to the neutron one, shown in 
Figs. 4 and 6. This does not occur for the d flavor, 
and the dominance of the neutron contribution are not 
spoiled by increasing A 2 as for the u flavor (cf. Fig. 6). 
This happens also because one half of the d content of 
3 He comes from the neutron, while only one fifth of the 
u one comes from it. 

Of course, the shape of the curves shown for GP M is very 
dependent on the nucleonic model used as input in the 
calculation, i.e., in this case, that of Ref. [25|. One should 
not forget that the aim of this analysis, for the moment 
being, is that of getting a clear estimate of the proton 
and neutron contributions to the nuclear observable, a 
feature rather independent on the nucleonic model. To 
demonstrate this property, we have plotted, in Fig. 7, 
the ratio of the proton to neutron contribution to & M \ 



r^ n (x, A 2 ,0 = &m(x, A 2 , / G m (x, A 2 , , (23) 

calculated using as input the G p ^' q corresponding to the 
model of Ref. [25| , and those of the very different model 
of Ref. [29| . The two ratio are slowly varying as a func- 
tion of x and are very close to each other, demonstrating 
the very weak model dependence of this feature of the 
result. 

Our IA calculation, so far, shows that, at very low 
—A 2 ~ 0.1 GeV 2 , there is a clear dominance of the 
neutron contribution on the proton one, and that such 
a dominance, even stronger for the d flavor, does not 
depend on the nucleonic model used in the calculation. 
Anyway, a couple of items have still to be investigated: 
i) even if the neutron contribution is dominating, the ex- 
traction of the neutron GPDs from it may be nontrivial 
and a proper strategy has to be studied; ii) although the 
very low A 2 values define the most interesting region, 
where, for example, the JSR can be checked, this region 
may be difficult to be reached experimentally. A proce- 
dure able to take into account the nuclear effects arising 
in an IA description also at higher A 2 would be very 
helpful. 

A positive answer to these two remaining problems will 
be given in the next section. 



EXTRACTING THE NEUTRON INFORMATION 
FROM 3 HE DATA 

Let us discuss the most relevant issue. 
It is convenient to rewrite our main equation, Eq. {T 
in a different form. 



A 2 
[GeV 2 ] 


s-i3,p, point 

Avl8 


si3,p, point 

Avl4 


^3,72, point 

Avl8 


s-i3,n, point 

Avl4 





-0.044 


-0.049 


0.879 


0.874 


-0.1 


0.040 


0.038 


0.305 


0.297 


-0.2 


0.036 


0.035 


0.125 


0.119 



TABLE I: The proton and neutron contributions to the mag- 
netic point like ff of 3 He, corresponding to the Avl8 [2j| and 
Avl4 24] nuclear interactions, for three low values of —A 2 . 



Let us define a variable z as follows: 

_ Ma t = _ Ma <- 
M £' Z M 



(1 + 0-3 



M A p+ 



M AP + 

where use has been made of Eqs. (|5J) - 
can be written therefore in the form 

M A ,M dz 



(24) 

(Hill). Eq. (H3D 



G^(x,A 2 ,0 = £ / - 9 %(z,A 2 ,0 



N Jx 3 
qN,q ( X a2 £ 



z z 



(25) 



where the off-diagonal spin-dependent light cone momen- 
tum distribution 



,0 . f iE f, 



dpP^(p,p + A) 



xt , c M AP + 

x 6 ^ z+ ^-ww 



(26) 



has been introduced and the quantity 

P%{p,p + A) = P?_, + _(p,P',E) 
- P?_,_ + (p,p',E) 



(27) 



for the sake of clarity, has been defined. 

Let us now perform the x— integral of the function 
G 3 ^(x, A 2 ,£)> given by Eq. ([25]). One obtains: 



/ 



dx G% q (x,A 2 ,0 = 

= e/^/t^' a2 '^ 9 (i' a2 '! 

= E / dx'G^(x',A 2 ,C) [ dzg%(z,A 2 ,0 
= Y,G N M q {^)Gl Wnt {A 2 ) 



N 

= Gl\A 2 ) . 



(28) 



In the equation above, G^?(A 2 ) is the contribution, of 
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FIG. 8: (Color online) The magnetic point like ff of 3 He (full), 
given by the sum of the neutron (dashed) and proton (dot- 
dashed) contribution (Eq. ([29]0. 

the quark of flavor q, to the nuclear ff; G^J q (A 2 ) is the 
contribution, of the quark of flavour q, to ff of the nucleon 
TV; G^ ,pmnt (A 2 ) is the so-called 3 He magnetic "point 
like ff", which would represent the contribution of the 
nucleon N to the magnetic ff of 3 He if N were point-like. 
The latter quantity has a relevant role in our discussion; 
we stress that it is given, in the present IA approach, by 

G 3 ! £ T ' point (A 2 ) = JdE J dpP%{p,p + K,E) 

= J dzg%{z,A 2 ,Z) . (29) 

This quantity, obtained in our Avl8 framework, is shown 
for N = p,n, and for their sum, G^ om£ (A 2 ) in Fig. 8. 
Obviously, this one-body property can be obtained from 
the wave function only and, at least for the low values 
of A 2 , it depends very weakly on the nuclear interac- 
tion. We have performed the calculation using also the 
Avl4 interaction, including the Coulomb repulsion be- 
tween the protons. The curves we have obtained cannot 
be distinguished from the Avl8 ones shown in Fig. 8, and 
have not been reported. The different magnetic point like 
ffs corresponding to the different nuclear interactions are 
therefore reported in Tab. 1 for three low values of —A 2 . 
Now it comes an important observation. The variable z, 
at the low values of A 2 and £ which are of interest here, 
is very similar to the light cone momentum fraction 
and the function g%{z, A 2 , £) very similar to a standard, 
forward light-cone momentum distribution. Besides, the 
nucleon dynamics in 3 He is, to a large extent, a NR one. 
This makes any light-cone momentum distribution, po- 
larized or unpolarized, strongly peaked around z = 1 
(see, i.e., the discussion in Ref. [12]). Evaluating Eq. 
(24) close to z = 1 means, as a matter of facts, that all 
the nuclear effects, but the ones due to the spin structure 
of the target, are negligible. Basically, this means that 



the momentum and energy distributions of the nucleons 
do not affect the result. If this is the case, from Eq. ([25]) 
one gets: 

N 

rM A /M 

x / dzg 3 N (z,A 2 ,0 + O(z-m) 

J X 

and, considering that, for 3 He, up to x ~ 1 one has 
g%(z, A 2 , £) ~ 0, being this function strongly peaked 
around z = 1, the lower integration limit in z can be 
put equal to for the x values relevant here (x < 0.7). 
In other words, at low values of A 2 and £, the following 
approximation of Eq. ([T3]) should hold: 

G 3 *(x,A 2 ,0 ~ ]TG^(z,A 2 ,£) 

N 

pMa/M 

x / dzg%{z,A 2 ,£) 
Jo 

= G 3 ™° int {A 2 )G™{x,A 2 ,® 

+ G 3 A f' pom *(A 2 )G^( a; ,A 2 ,C).(31) 

Clearly, in the last line, use has been made of Eq. ([29]) . 
If the factorized expression above were a good approxi- 
mation of Eq. ([13]) . it would be very helpful. In facts, 
Eq. ([3T]) is very simple: all the nuclear effects are hidden 
in the magnetic point like ffs, quantities which can be 
obtained directly from the wave function, without con- 
sidering the complicated spectral effects described by Eq. 
([13]) : moreover, and more important, these quantities, at 
low —A 2 , are theoretically well known and, as Tab. 1 
shows, depend very weakly on the used nuclear interac- 
tion. 

Figs. 9 and 10 demonstrate that Eq. ([3T]) approxi- 
mates nicely the full result Eq. ([T3]) in the kinematics of 
interest here, the two quantities differing of few per cents 
at most for x < 0.7. This is an important point: Eq. 
([3T]) can now be solved to extract the neutron GPD G 7 ^: 

G^ tr (x,A 2 ,0 * . * {G 3 M (x,A 2 ,0 

M \^ ) 

- G 3 ^ poM (& 2 )G p M (x,A 2 ,0}m) 

an equation which could be used to obtain the neutron 
information from 3 He and proton data, simply correcting 
by means of the well known magnetic point like ffs. 

We have checked the validity of the proposed extrac- 
tion procedure by evaluating Eq. ([32]) . using for the mag- 
netic point like ffs the ones corresponding the Avl8 in- 
teraction, for the proton GPD G P M the one given by the 
model [lij , for & M the one calculated by means of Eq. 
([T3]h i.e., we are simulating 3 He data by using our best 
calculation. If the extraction procedure were able to de- 
scribe exactly, in an effective way, the nuclear corrections 
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FIG. 9: (Color online) The quantity x 3 G 3 M (x, A 2 ,£) of 3 He 
in the forward limit, together with the proton and neutron 
contribution (full lines), compared with the approximations 
to these quantities given by the factorized form Eq. (|3T|) 
(dashed lines). 



FIG. 11: (Color online) The quantity x 3 G^ q (x, A 2 , f ) for the 
neutron in the forward limit, evaluated in the model of Ref. 
[25| and labelled n, given by the sum of the u, labelled n u , 
and d, labelled rid contribution (full lines), compared with the 
approximations to these quantities given by Eq. (|32|) (dashed 
lines) (see text). 



-0.2 




neutron 



x 3 G 3 M (x,A 2 ,£) 



0.2 



0.4 



0.6 



0.8 

x 3 



FIG. 10: (Color online) The same as in Fig. 9, but at A 2 
-0.1 GeV 2 and & = 0.1. 




FIG. 12: (Color online) The same as in Fig. 11, but at A 2 
-0.1 GeV 2 and f 3 = 0.1. 



predicted in IA, the obtained G^f x r should be equal to 
the neutron quantity, , used as input in the calcula- 
tion, i.e., again, the one predicted for the neutron by the 
model of Ref. U^. 

Figs. 11 and 12 demonstrate that these two quanti- 
ties, and the u and d flavor contribution to them, differ 
at most by few percents in all the relevant kinematical 
range. This means that the growth of the proton con- 
tribution to the 3 He observable, which seemed to hinder 
the extraction of the neutron information, in particular 
for the u flavor, is governed by the —A 2 behavior of the 
magnetic point like ffs, quantities which are under the- 
oretical control. By using them, the extraction of the 
neutron GPDs, close to the forward limit, is safe and 
basically model independent. 

These feature is even more evident in Fig. 13, where 



the ratio 

r n (x,A 2 ,0=G n M extr (x,A 2 ,0/G n M (x,A 2 ,0 , (33) 

is shown. Obviously, this ratio would be one in case the 
extraction procedure were working perfectly. It is clear 
that, in all the relevant kinematics, the extraction proce- 
dure takes into account the nuclear effects introduced in 
the I A analysis up to x ~ 0.7 with an accuracy of a few 
percents. 

Besides, the quality of the extraction does not depend 
on the nucleonic model used. In Fig. 14, in one particular 
kinematics, it is shown the ratio r n , Eq. ([33]) . obtained 
using the two different models of GPDs considered in this 
work, i.e. the ones ofRefs. jHHil. It is clearly seen that 
the procedure works by using both models within a few 
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FIG. 13: (Color online) The ratio Eq. ((33), in the forward 
limit (full), at A 2 = 0.1 GeV 2 and £3 = (dashed) and at A 2 
= 0.1 GeV 2 and £3 = 0.1 (dot-dashed). 
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FIG. 15: (Color online) The same as in Fig. 3, but at A 2 
=-0.2 GeV 2 and & = 0.2. 




FIG. 14: (Color online) The ratio Eq^ (J33J) , at A 2 0.1 GeV 2 
and £3 = 0, using the model of Ref. [25| for the nucleon GPDs 
(dashed) and the one of Ref. [29[ (full) to evaluate Eq. (|13p . 



percent accuracy. 

Let us show now that the technique works at higher 
—A 2 values (but one should keep in mind that, for the 
procedure to be reliable, the IA should work). Figs. 15- 
18 are devoted to show this. All the steps are repeated, 
in a sort of useful summary of the procedure. In Fig. 15, 
the proton and neutron contributions to 3 He are shown 
at A 2 =-0.2 GeV 2 and £3 = 0.2, and they are found to 
be comparable in size; in Fig. 16, it is shown that this 
problem is serious, in particular, for the u flavor. In Fig. 
17 and 18 the good quality of the extraction procedure 
is demonstrated even in this less forward situation, by 
showing the extracted neutron GPD and its ratio to the 
model used as input in the calculation, respectively. This 
is a good piece of news for the experimental programme, 
in case that extremely small values of —A 2 could not 
be reached. However, we reiterate that the very good 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 Oi 



x 3 

FIG. 16: The same as in Fig. 6, but at A 2 =-0.2 GeV 2 and 
6 = 0.2. 



extraction scheme is really useful where IA provides a 
reasonable description of the process, a region which co- 
incides with -A 2 < 0.15 GeV 2 . 

In closing this section, let us mention that another pro- 
cess to access the neutron GPDs is incoherent DVCS off 
the neutron in nuclear targets, i.e., the process with the 
struck neutron detected in the final state, in coincidence 
with the scattered electron and the produced photon. An 
experiment of this type will be performed at the 12 GeV 
program of JLab [33] for a 2 H target. We plan therefore 
to investigate also incoherent DVCS off the neutron in 
3 He, although these kind of processes could be spoiled 
by Final State Interactions of the detected neutron. 



CONCLUSIONS 

In this work we have thoroughly described an IA cal- 
culation of the GPDs H q and E q of 3 He. In proper limits, 
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FIG. 17: (Color online) The same as in Fig. 11, but at A 2 
=-0.2 GeV 2 and f 3 = 0.2. 
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the correct constraints are recovered. Coherent DVCS off 
3 He at low A 2 turns out to be strongly dominated by the 
neutron contribution, in particular for the d flavor. From 
such a neutron contribution, the neutron information has 
to be extracted. A procedure has been described to take 
into account the nuclear effects included in the IA anal- 
ysis and to safely disentangle the neutron information 
from them, even at moderate values of the momentum 
transfer. The only theoretical nuclear ingredients are the 
neutron and proton contributions to the magnetic point 
like form factors of 3 He, quantities which are under good 
theoretical control at low momentum transfer and that 
encode correctly all the nuclear effects described in an 
I A framework. If high values of A 2 were reached, the I A 
description would not be reliable and two-body currents 
and three-body forces would have to be included into the 
approach. We have checked that the dominance of the 
neutron contribution and the proposed extraction proce- 
dure do not depend on the model of the nucleon GPD 
used in the calculation. Our results confirm strongly co- 
herent DVCS off 3 He at low momentum transfer as a key 
experiment to access the neutron GPDs. It will be very 
interesting and useful to perform a Light-Front analysis 
of the process, which already started in SiDIS (Hf, so to 
have, from the beginning, a relativistic framework for the 
investigation. 
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APPENDIX A: GPDS FROM THE COMPONENTS OF THE LIGHT CONE CORRELATOR 

In this work, both the nucleus (because the relevant momentum transfer is low) and the nucleon (whose dynamics is 
governed by the Schrodinger equation) are treated non-relativistically (NR). A NR expression of the spinors appearing 
in Eq. (pQ) has therefore to be used, in order to find explicit relations between the correlator F^f ,fl and the GPDs. 

In the NR limit one has (see, e.g., [HI): 

/ 



u(p, s) 



and: 



Vp^Xs I NR 



2m 



Xs 



2rn 



Xs 



(34) 



u(p \s f )rfu(p,s) 



Zrnxlxs , 

2m X \, i^^ k ) Xs . 



(35) 
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Properly choosing the components of 7 M , it is possible to find the most convenient relation between F^f'^ and the 
GPDs. By taking for \ the eigenstates of a 3 , being z the quantization axis, fixing \i — 0, one obtains: 

F^Or.A 2 ^) = -^H*{x,A\0 

?Q,A,0, A 2 ^ _ P ?,A,0/ A 2 



Fl^(x, A 2 ,£) = F^(x,A\0 = . (36) 



In the same way, for \i = 1 one gets 



F^\x,A\i) = 0, 



+ + 

A3 

*Z,A,1/ A 2 >n _ _ rpq,A,l/ A 2 >x _ A 2 



F^\x,^,i) = -F^(x,A^,e) = -^+G^(x,A^,e). (37) 



2p 

Collecting all the above results, one gets: 



H?(x,A*,t) = —Fl'f' (x,A 2 ,£) 



A 3 

which coincides with Eq. (fTTj) . 



G^(x,A 2 ,0 = ^F^(x,A 2 ,0, (38) 



APPENDIX B: THE LIGHT CONE CORRELATOR FROM NUCLEAR OVERLAPS IN IA 



The overlap Eq. ([T5]) is the crucial nuclear ingredient to evaluate the light cone correlator, Eq.(l), for 3 He. Let us 
show how it is obtained through the quantity at our disposal, which is actually the following: 



I^iyuE) = /"d% {XM x ,L p M p \ l -J z )Yl pMp ^l ill ) 

J M p ,Mx 

x E0'i^i2,^I^M x )( yi) ^;^ mi2 ' s |*y;). (39) 

Here above, p denotes quantum numbers of the struck nucleon, treated as a free particle in IA, with spin projection a, 
isospin projection r, orbital angular momentum and its third component L p , M p , respectively. '0^ j 2 r J mi2S represents 
the two-body recoiling system with total angular momentum and its third component ji2, mi2, spin £ and excitation 
energy E. rj is an additional quantum number necessary to describe the OAM mixing in the interacting two-body state 
in the continuum. ji2 is coupled to the spin | of the interacting nucleon to give the intermediate angular momentum 
X. The description of the three-body system is based on the Avl8 calculation of the wave function of Ref. [22], given 
in terms of the following Jacobi coordinates yi and xi: 




yi 

xi 



2 

"7! 

[r 2 - r 3 



ri 



r 2 + r 3 



(40) 



where (with i = 1,2, 3) are the CM variables of the nucleons in the trinucleon, constrained by i*i + r2 + rs = 0. The 

I j 

trinucleon bound state, ^ 3 2 ^(xi,yi), can be written in terms of the variables yi and xi using the following basis for 
the isospin-spin- angular part: 
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M p M x mi2(J 

X ^i 2 mi 2 (^i) (42^12,2^1^12^12) 



mi 2 /x 

X 



E<r,T 12 r 12 |ii) 77 T^l^M^Jxf (41) 

through the overlaps 



\ 2 , ^'22' 2 

rri2 



ii^dxil, | yi |) = y dfi fil y [^^(ii.yi)] *?£(xi,yi) (42) 
with the constraints 

(_l)«/+ip = + i 

(_l)</+s+Ti 2 = _! ; (43) 

following from the parity of the nucleus and the Pauli principle, respectively. 

The quantity at our disposal, Eq. ([39]) . can be defined in terms of the overlaps Eq. ([43]) . as follows: 

iL^iyuE) = <ir,r 12 r 12 |ii> I ^"(lyil*^ (44) 

where u J g™ ,E,ri \x\) is the two-nucleon radial function. 

Using this formalism to calculate the spectral function, Eq. (Q , one has to obtain the product of two intrinsic 
overlaps, integrated over the angles of the relative motion of the recoiling pair. This quantity is expressed through 
the overlaps Eq. ([44]) which, as previously said, are the quantities at disposal, as follows: 



t i 



J dn t (P'S'\~ts u f's f )(ps,ts t \PS) = 



= [ dp'e-**? [ dpe& E E EE n,M^p)^M'0V) 
x (XM x ,L,M\±J z )(j 12 m 12 ,±s\XM x )(J' z ±\L'M',X'M x ,) 

x (X'M x ^sU 2 m 12 )jf;*' M ^ (45) 

Using these relations, together with Clebsch-Gordan coefficients and spherical harmonics properties, the final expres- 
sions for the relevant components of the light cone correlator are found: 



x | dpp 2 j Lp (pp)if^ j ^(p,E) 

/l /.27T 
dcosO p J d(j) p PL p (cos9^ 



J dp'p'*j Lp (\p + A\p'^^ . (46) 
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x i v (|p + A|p')jL p (w)V(^^)// 3 (p,t)r Lp M p (^)(-i) i "' + 
x (F e r L * p/>Mp _ 2 (fi p+A ) - £> e y L * p>> M p (n P +A))F^ 1 ( a /^', a 2 ) 



L.+L p 



(47) 



where 7' = L p , M p , X, j'12, S, A, L p /, X 7 . and are terms containing the product of four Clebsch-Gordan coeffi- 
cients, depending on the set of quantum numbers g = L p >L pi M p , X, X 7 , J12. 

Using the latter equations in Eqs. JTTJ, Eqs. (fT2j) and (fT3|) can be eventually evaluated, through the overlaps at 
our disposal. 

(1993). 

[17] M. Diehl, Phys. Rept. 388, 41 (2003). 
[18] A.V. Belitsky and A.V. Radyushkin, Phys. Rept. 418, 1 
(2005). 

[19] S. Boffi and B. Pasquini, Riv. Nuovo Cim. 30, 387 (2007). 
[20] L. L. Frankfurt and M. I. Strikman, Phys. Rept. 76 

(1981) 215; Phys. Rept. 160 (1988) 235. 
[21] A. Kievsky, E. Pace, G. Salme and M. Viviani, Phys. 

Rev. C 56, 64 (1997). 
[22] A. Kievsky, M. Viviani and S. Rosati, Nucl. Phys. A 577 
(1994) 511. 

[23] R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. 

Rev. C 51 (1995) 38. 
[24] R. B. Wiringa, R. A. Smith and T. L. Ainsworth, Phys. 

Rev. C 29, 1207 (1984). 
[25] I. V. Musatov and A. V. Radyushkin, Phys. Rev. D 61 
(2000) 074027. 

[26] A. V. Ra dyushkin, Phys. Lett. B 449, 81 (1999) 

pep-ph/9810466] . 
[27] M. Gari and W. Krumpelmann, Phys. Lett. B 173, 10 
(1986). 

[28] M. Diehl, Th. Feldmann, R. Jakob and P. Kroll, Eur. 
Phys. J. C 39, 1 (2005); M. Guidal, M. V. Polyakov, 
A. V. Radyushkin and M. Vanderhaeghen, Phys. Rev. D 
72, 054013 (2005). 
[29] S. Scopetta an d V. Vento, Eur. Phys. J. A 16, 527 (2003) 

[hep-ph/0201265] , 
[30] L. E. Marcucci, D. O. Riska and R. Schiavilla, Phys. Rev. 

C 58, 3069 (1998). 
[31] F. Baroncini, A. Kievsky, E. Pace and G. Salme, AIP 

Conf. Proc. 1056 (2008) 272. 
[32] A. Amroun, V. Breton, J. M. Cavedon, B. Frois, 
D. Goutte, F. P. Juster, P. Leconte and J. Martino et 
a/., Nucl. Phys. A 579, 596 (1994). 
[33] JLab Exp. E1211003, Spokesperson: S. Niccolai. 
[34] S. Scopetta, A. Del Dotto, E. Pace, G. Salm e, II Nuovo 

Cimento C, 35, 101 (2012) larXiv: 11 10.67571 [nucl-th]. 
[35] M.E. Peskin, D.V. Schroeder, "An introduction to Quan- 
tum Field Theory", Addison- Wesley Advanced Book 
Program, New York, 1994. 



[1] D. Mueller e t al, Fortsch. Phys. 42, 101 (1994) 

[hep-ph/9812448] . 
[2] A. V. Radyushkin, Phys. Lett. B 380, 417 (1996); Phys. 

Rev. D 56, 5524 (1997). 
[3] X. -D. Ji, Phys. Rev. Lett. 78, 610 (1997); Phys. Rev. D 

55, 7114 (1997). 
[4] M. Burkardt, Phys. Rev. D 62, 071503 (2000) [Erratum- 

ibid. D 66, 119903 (2002)]. 
[5] M. Vanderhaeghen, P. A. M. Guichon and M. Guidal, 

Phys. Rev. D 60 (1999) 094017. 
[6] HERMES Collaboration, A. Airapetian et al, Nucl. 

Phys. B 829, 1 (2010); Phys. Rev. C 81, 035202 (2010); 

M. Mazouz et al, (JLab Hall A), Phys. Rev. Lett. 99, 

242501 (2007). 

[7] M. Guidal, Phys. Lett. B 693, 17 (2010); Phys. Lett. B 

689, 156 (2010). 
[8] E. R. Berger, et al. Phys. Rev. Lett. 87 (2001) 142302; 

F. Cano and B. Pire, Eur. Phys. J. A 19 (2004) 423. 
[9] A. Kirchner and D. Mueller, Eur. Phys. J. C 32 (2003) 

347; 

[10] M. V. Polyakov, Phys. Lett. B 555 (2003) 57; V. Guzey, 
M. Siddikov, J. Phys.G32, 251 (2006); V. Guzey and 
M. Strikman, Phys. Rev. C 68, 015204 (2003); S. Liuti 
and S. K. Taneja, Phys. Rev. D 70 (2004) 074019; Phys. 
Rev. C 72 (2005) 032201; Phys. Rev. C 72 (2005) 034902; 
V. Guzey, Phys. Rev. C 78, 025211 (2008); V. Guzey, 
A. W. Thomas and K. Tsushima, Phys. Lett. B 673 
(2009) 9. 

[11] J. J. Aubert et al. [European Muon Collaboration], Phys. 

Lett. B 123, 275 (1983). 
[12] S. Scopetta, Phys. Rev. C 70 (2004) 015205; 

Phys. Rev. C 79 (2009) 025207. 
[13] M. Rinaldi and S. Scopetta, Phys. Rev. C 85, 062201(R) 

(2012). 

[14] J. L. Friar, B. F. Gibson, G. L. Payne, A. M. Bernstein 
and T. E. Chupp, Phys. Rev. C 42 (1990) 2310. 

[15] C. Ciofi degli Atti, S. Scopetta, E. Pace and G. Salme, 
Phys. Rev. C 48 (1993) 968. 

[16] R. -W. Schulze and P. U. Sauer, Phys. Rev. C 48, 38 



